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We review the current status of theoretical predictions for mixing quantities and lifetimes in 
the B-sector. In particular, due to the first non-zero measurement of the decay rate difference 
in the neutral -B s -system, AF S /F S = 17.6%±2.9% by the LHCb collaboration and very precise 
data for tb b from TeVatron and LHCb our theoretical tools can now be rigorously tested and 
it turns out that the Heavy Quark Expansion works in the B-system to an accuracy of at least 
30% for quantities like Fi2, which is most sensitive to hypothetical violations of quark hadron 
duality. This issue that gave rise in the past to numerous theoretical papers, has now been 
settled experimentally. Further data will even allow to shrink this bound. For total inclusive 
quantities like lifetimes the compliance is even more astonishing: r^f Cb /rfJ AG = 1.001±0.014 
is in perfect agreement with the theory expectation of rs s /r_B d = 0.996. ..1.000. 
Despite the fact that the new data show no deviations from the standard model expectations, 
there is still some sizable room for new physics effects. Model-independent search strategies 
for these effects are presented with an emphasis on the interconnection with many different 
observables that have to be taken into account. In that respect a special emphasis is given to 
the large value of the di-muon asymmetry measured by the DO collaboration. 



1 Introduction 

RecentljQ several hints for new physics effects in the flavour sector, see e.g. I ^ I ^H^ triggered a 
lot of interest. In particular there were several independent experimental indications^SIEIZIHMIO] 
for a large CP violating phase in ^-mixing, wh ich would be in clear contradiction to the tiny 
standard model (SM) expectation ^ 1 | 12 | 13 | 14 | 1 5] 

Unfortunately the huge expectations in spectacular new physics effects were not c onfirmed by 
precise measurements at LHCb or by follow-up measurements at the TeVatron, e.g. | 16 | 17 | 18 119| 
The current disappointment about this situation can be nicely visualized with Fig.JTJ). § 



"Before the Lepton Photon Conference in August 2011. 

b A r ece nt review on the implications of this new experimental results on several new physics models can be 
found in^4 




Figure 1: How the new LHCb results might appear at first sight. Previous measurements at TeVatron raised big 
expectations in large new physics effects, which were unfortunately not confirmed by more precise measurements 
at LHCb. It is, however, important to note that we gained many theoretical insights by the first measurement of 
AF S and by other very precise measurements. This will be elaborated in this paper. 



One of the main aims of this talk is to emphasize the fact that we now have for the first time 
a measurement of a non-zero value of the decay rate difference in the neutral -B s -system, ATp^l 
This quantity raised a lot of interest in the last thirty-five years and several open questions 
related to the theoretical determination of AT S have now been answered experimentally. Another 
important point is the fact that there is still some ample space for new physics effects in B- 
mixing. To distinguish possible new physics effects from hadronic uncertainties, however, a 
considerably higher theoretical precision for questions like penguin pollution is now necessary, 
compared to what would have been needed in the case of a huge mixing phase. 
In Section 2 we give a short introduction to the mixing formalism and to its experimental 
evidence. We also briefly touch the calculational tools, in particular the Heavy Quark Expansion 
(HQE). In Section 3 historic attacks to the HQE and their experimental resolution are described. 
Therefore we discuss theoretical arguments in favour of the HQE as well as the missing charm 
puzzle, the lifetimes of the A^-baryon and the i? s -meson, the decay rate differences ATj and Ar s 
and semi leptonic CP asymmetries. This chapter contains the main results: the proof of the 
applicability of the HQE to determine T\2- Model independent investigations of possible new 
physics effects in mixing are presented in Section 4, while Section 5 concentrates on new physics 
contributions and its constraints to the absorptive part of the mixing diagram, T^. Such effects 
were discussed in order to explain the large central value of the di-muon asymmetry measured 
by the DO collaboration. Since it turned out experimentally that new physics effects in S-mixing 
are not huge, investigations of the remaining space for effects beyond the standard model require 



now a much better theoretical control of e.g. penguin contributions, this is discussed in Section 
6. Finally we conclude. 



2 Mixing Formalism 



2. 1 Basic mixing formulae 

The time evolution of a decaying particle B{t) = exp [— imst — Ts/Zt] with mass and 
lifetime tb = 1/Tb can be written as a differential equation. In the case of a two-state system 
(e.g. the neutral meson B q = (bq) and its anti-particle B q = (bq)) one expects a diagonal mass 
matrix M q and a diagonal decay rate matrix T q : 
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It is known since a long time (e.g. that due to weak interactions transitions like B^ s — > 
Bd,s are possible in the neutral 5-system. They are triggered by the so-called box diagrams. 
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The box diagrams lead to off-diagonal elements in the mass matrix M q (denoted by M\ 2 
to off-diagonal elements in the decay rate matrix T q (denoted by T^ 2 ). Mf 2 corresponds to the 
part of the box-diagram, where the internal quarks are off-shell (e.g. the internal top-quark), 
while r^ 2 denotes the part where the internal quarks are on-shell (in the diagrams shown above 
this can only happen for the up- and charm-quark). Both Mf 2 and T\ 2 can be complex due to 
their dependence on CKM-elements. 

Diagonalization of M q and T q gives mass eigenstates B q ^ (H=heavy) and B q ^ (L=light) with 
the masses Mi, M? and the decay rates T g H , Tf , see e.gP^i 



B q ,H ■= V B q - q B q 



B, 



q ,L 



pBg + qBq with \p\ 2 + \q\ 2 



(2) 



In the limit of no CP violation in mixing, the heavy mass eigenstate is CP-odd and the light one 
CP-even. Theoretical calculations of the box diagrams (Mf 2 and T^ 2 ) can be related to three 



12/) 



observables, expressed in terms of |M-f 2 |, |r^ 2 | and the relative phase 4> q = arg(— M-^/T^ 
which describes CP-violation in mixing; the individual phases of Mf 2 and T^ 2 are unphysical: 

• Mass difference: 



The mass difference between the heavy and the light mass eigenstates of the neutral B q 
meson is given by the modulus of the dispersive part of the box diagram. 

AM, := Mjj-Ml = 2|M? 2 | - sin 2 ^ + ...j . (3) 

In the standard model | r ^ 2 1 / 1 JViT ^2 1 is °f the order of 5 per millet™ for the .Bg-system. 
Therefore the corrections in Eq.([3|) are completely negligible. Since M\ 2 describes the 
virtual part of the box diagram, it is sensitive to heavy internal particles like the top- 
quark, but also to hypothetical new states like SUSY-particles, KK-states, ... . So, the 
mass difference is expected to be very sensitive to new physics effects. 



Decay rate difference: 



The decay rate difference between the heavy and the light mass eigenstates of the neutral 
B q meson is given by the modulus of the absorptive part of the box diagram. 

Ar g: =rl-r^ = 2|r? 2 |co S( A^i + i^^sin 2 ^ + ..^ . (4) 

In the standard model <j> s is close to zer so the cosine gives to a good approximation 

a value of one. Since T\ 2 describes the on-shell part of the box diagram it is sensitive to 
light (with a total mass below Mb) internal particles, like the up- and charm-quark. As 
we do not expect sizable new physics contributions to CKM-favoured tree-level decays like 
b —> ccs, it seems reasonable to assume almost no (i.e. below the hadronic uncertainties) 
new physics effects in Tf 2 . Therefore the only way for new physics to affect Ar s is the 
phase <p s . In the standard model 4> s is close to zero, so new physics can only lower Ar s 
compared to its standard model value^l. We will challenge, however, the assumption that 
Y\ 2 can not be affected by new physics in Section 5. 

Flavour specific /semi leptonic CP asymmetries: 



A flavour specific decay B q — > f is defined by the property 

— The decays B q — > / and B q — > / are forbidden. 

To simplify the analysis of untagged decays one often demands in addition that no direct 
CP violation occurs in the decay, i.e. \(f\B q )\ = \(f\B q )\. Examples of such a decay are 
B s — > D~ir + or the semi leptonic decay B q — > Xlv - hence also the name semi leptonic 
CP asymmetry. The asymmetry is defined as 





a sl = a fs 
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Mf 2 



T(B q (t)^f)-T(B q (t)^f) = _ 2 ( g 
T(B q (t)^ f) + T(B q (t)^ /) V P 

T q AT 

sin^ = Im^ = £^tanV (5) 



In the standard model these asymmetries are suppressed by the small values of |r^ 2 |/|M^ 2 | 
and 4> q . Measuring a sizable value for the semi leptonic CP asymmetries would be a clear 
indication for new physics, provided that we have a precise control over the theoretical 
framework. Theoretical arguments in favour of this control and the experimental proof of 
this control, which is now available, will be discussed in detail below. 

The phenomenon of particle- antiparticle mixing is a macroscopic quantum effect, which is by 
now well established in several systems of neutral mesons. 

1955 X°-system: Mixing in the neutral if-system was theoretically developed in 1955 by Gell- 
Mann and Pais^S. Based on that framework the phenomenon of regeneration was predicted 
in the same year by Pais and Piccionil^]. Experimentally regeneration was confirmed in 
1960 A huge lifetime difference between the two neutral K- mesons was established 
already in 1956^^21. 

1986 ^-system: Mixing in the S^-system was found 1986 by UA1 at CERN^ (uAl attributed 
the result however to B* mixing) and 1987 by ARGUS at DES¥^. The large result for the 
mass difference AM^ can be seen as the first clear hint for an (at that time) unexpected 



large value of the top quark mas;P^ |. For the decay rate difference currently only upper 
bounds are available, see^ZI f or the most recent and most precise bound. 

2006/12 i^-system: The large mass difference in the -B s -system was established by the CDF col- 
laboration at TeVatron^3. In 2012 the LHCb Collaboration measured for the first time a 
non- vanishing value of the decay rate difference in the £? s -system^H 

2007 Z?°-system: Here we have several experimental evidences (BaBar, Belle, Cleo, CDF, E791, 
E831) for values of AT/T and AM/T at the per cent level, but we still do not have a 
single measurement with a statistical significance of more than five standard deviations. 
The combination of the datJ3 shows, however, unambiguously that D mesons mix. 

A more detailed discussion of the experimental evidence for mixing can be found e.g. in^Sl. 
Because of several reasons, mixing observables seem to be very well suited for the search of 
hypothetical new physics effects. First, these observables are experimentally established - as 
summarized above. Next, in the standard model the box diagrams are strongly suppressed, 
because of being fourth order in the weak interaction. So even small new physics effects might 
give comparable and therefore detectable contributions. Finally, hadronic effects are well under 
control in the B mixing observables. This will be elaborated below. 



2. 2 Theoretical determination of M\2 

Calculating the virtual part of the box diagrams with internal up-, charm- and top-quark and 
using the unitarity of the 3x3 CKM matrix one finds that the by far dominant contribution is 
given by the top-quark. It reads 

M q 12 = ^{V t * q V th ) 2 M^S Q {x t )B Bq f 2 B M Bq f, B . (6) 

The evaluation of the 1-loop box-diagram gives the so-called Inami-Lim function 5p(xt)^^ with 
xt = m 2 /Myy- Perturbative QCD corrections to the box-diagrams are denoted by they 
turned out to be ample. In the end one also has to determine the non-perturbative matrix 
element of a four quark operator Q (the four external legs of the box diagram) switched in 
between the states of the physical B and B mesons. For historical reasons this matrix element 
is parameterized in terms of a bag parameter B Bq and a decay constant f Bq : 

(B g \Q\B g ) = S -B Bq f 2 B M 2 Bq . (7) 
Q = ^7/,(l-75)&a-9/37 M (l- 75)^/3 (8) 

a and j3 denote colour indices. The bag parameter and the decay constant have to be determined 
with non-perturbative methods like lattice QCD or QCD sum rules. 

Formally one has performed in Eq.([6]) an operator product expansion (OPE), by integrating 
out the heavy VF-boson and the heavy top-quark. Corrections to that expansion are expected 
to be of the order m 2 /M^, i.e. smaller than four per mille and thus far below the hadronic 
uncertainties. 

Comparing experiment (HFACp^) with theory^^ one sees a very nice agreement, but also the 
fact that the theoretical error is considerably larger than the experimental one. The theory error 



c To avoid a very large value of the top quark mass, also different new physics scenarios were investigated, in 
particu lar a scenario with a heavy fourth generation of fermions and a top quark mass of the order of 50 GeV, 
see e.g.EH. 



is currently dominated by the hadronic parameters. 



AA/| M = 0.543 ± 0.091 ps _1 , 



AM, 



SM 



17.30 ± 2.6 ps 



-i 



AM d = 0.507 ± 0.004 ps" 1 
(ALEPH, CDF, DO, DELPHI, L3, 
OPAL, BABAR, BELLE, ARGUS, CLEO), 
AM 8 = 17.69 ± 0.08 ps" 
(CDF, DO, LHCb). 



0) 



(10) 



For the theory prediction in Eq. © and Eq. (|10p taken fronjm a very conservative average^ of 
lattice determinations for the matrix element in Eq. ([7]) was taken as an input 0. In particular 



f Ba = 231 ± 15 MeV 



(11) 



was used for the d ecay constant. Recently several new evaluations of the decay constant were 
performe d^ | 41 l 4 2] thi s new numbers one obtains: 



Reference 


fB s 


AM S 


HPQCD (1110.4510) 
Fermilab/MILC (1112.3051) 
CKMfitter new average (1203.0238) 


225 ± 4 MeV 
242 ± 10 MeV 
229 ± 6 MeV 


16.4 ± 1.0 ps" 1 

19.0 ± 1.8 ps" 1 

17.1 ± 1.7 ps" 1 



(12) 



For the first two results for AM S in Eq. (|12}l only the value and the error of the decay constant was 
changed compared to Eq. (|10p . For the third number also a new average for the bag parameter 
was taken into account. In view of th e relat ive large difference of the central values of the 
decay constants (17 MeV) determined prefer to stay currently with our conservative 

prediction given in Eq. (]10p . 

Here a further reduction of the theoretical errors and a future reduction of the difference between 
distinct evaluations will be very desirable, because the mass differences give important bounds 
on 

fitsSa 

of the unitarity triangle and also to searches for new physics effects. 



2.3 Theoretical determination of T 12 

The theoretical determination of rf 2 is much more complicated than the one of Mf 2 , sketched in 
the previous subsection. T\ 2 is sensitive to real intermediate states like the up- and the charm- 
quark. Thus one can not integrate out at once all particles inside the loop of the box diagram, 
one has to follow instead a two step procedure: 

1. OPE I: Integrate out the heavy W boson. This is similar to the case of M^ 2 . Again, 
corrections to that approximation are expected to be at most at the per mille level. 

2. OPE II: Peform an expansion in inver se powers of the heavy quark mass, the heavy quark 
expansion (HQE)l 44 l 4 5 l 4 « l 4 7 l « EQl. 

Then the expression for T\ 2 can be written in the following form: 



1 12 



AY 



A is a mass scale, which is expected to be of the order of Aqcd- Its precise value has, however, to 
be determined by explicit calculation. As will be shown below, for the 5 s -system the expansion 



d For a comparison of different wa ys to a verage the lattice results, see e.g. the discussion in^ 
Tor more recent averages, see e.gJ^^. 



parameter A/m^ is a little smaller than 1/5. The 1^ s are products of perturbatively calculable 
Wilson coefficients and non-perturbative matrix elements of operators with a mass dimension 
equal or larger than six. Such a matrix element is also parameterized in terms of the decay 
constant fs and a bag parameter Bk, i.e. 1^ oc /j| C-kBk- 

First estimates of T^ 0) were made starting from 1977 on, see e.g. I 51 | 52 | 53 | 54 | 55 | 56 | At that 
time no effective Hamiltonian and no non-perturbative determinations of the decay constant 
and the bag parameters were available. Since 1977 clearly the determination of T^- 1 improved 
a lot. The subleading l/m^-corrections were determined for the first time in 1996^3 for 
B s mesons, they turned out to be quite large. The NLO-QCD corrections for the leading 
CKM-structure in AT S were calculated already in 1998^1 In 20 03 this result was extended to 
the subleading CKM-structure by two independent groups^HH. These subleading terms give, 
however, the dominant contribution to the semi leptonic CP asymmetries and Ar^, while they 
are a small, but non-negligible contribution to the decay rate difference Ar s . In 2006^^ several 
theoretical improvements were incorporated in the theoretical determination of Tf 2 , like an use 
of different quark mass schemes or summing up large logarithms of the form m^/ml In w? c jm^ 
to all orders (following^!) as we n as a change of the operator basis. Calculating Ti2 to the order 
A 4 /m^ one gets 

Tt 2 = C(Q) q + Cs(Qs) q + C s (Qs) q + £ (c^Ri^ + C^iRi),) , (14) 

i=l,2,3 

with Wilson coefficients C, Cs, Cs and Cr and matrix elements of the operator Q appearing in 
M\2 (defined in Eq.(|8|)) and of two new operators 

Qs = 5a (1 + l5)b a ■ 5/9(1 + 75)&/3 , (15) 
Qs = 5a(l + 75)&/9 • 5>(1 + 75)&a • (16) 

Ri is an acronym for different operators of dimension seven, appearing in I^, which can be found 
e.g. m 1 ^ InE3it was noted that Q, Qs, Qs are n °t independent, they are related via 

R = Q S + ai Q s + \Q ■ (17) 

The coefficients have the form on = 1 + 0(a s (m&)) and Ro is of order l/m&. So, to leading order 
in 1/mb the quantity Y\ 2 can be expressed in terms of two independent operators, e.g. (Q, Qs), 
denoted as the old basis or (Q,Qs), denoted as the new basis. 

old basis: Tl 2 = C oM (Q) q + Ct d (Qs) q + C o ^(R ) q +J2 (c Ri (R i ) g + C Ri (R) q ) . (18) 

1=1,2,3 

new basis: T\ 2 = C°™{Q) q + C s cw (Q s ) q + Cf^(R ) q +£ (Cr^), + Cr^),) .(19) 

i=l,2,3 

In principle both expressions are absolutely equivalent. But, since the bag parameters of the 
different operators are not known with the same precision the choice of the basis has a sizable 
effect on the resulting accuracy in the determination of T\ 2 . The best known matrix element 
is the one of the operator Q. This operator also arises in the mass difference AM q and there 
are many non-perturbative determinations on the market. For the matrix element of Qs or Q s 
much less calculations are available, while for the matrix element of Rq only vacuum insertion 
approximation is used - this corresponds to an uncertainty of 50%! Looking now at the sizes of 
the coefficients in Eq. (fT8|) and Eq. ([T9|) one findlPl 



• ^iio™ ^ ^Ro : ^ n ^ e new basis the coefficient of the operator Rq is considerably smaller 
than in the old basis - this favours the new basis, because the matrix element of i?o has a 
very large uncertainty. 

• In the old basis Cg ld is the dominant coefficient, in the new basis the coefficient C new of 
the best known operator is dominant - this favours the new basis. 

• In the old basis, there is an accidental strong numerical cancellation in the coefficient C old , 
which can lead to an overestimate of the theoretical uncertainty in the old basis. 

• In the new basis the by far dominant contribution to the theory prediction for AT q / AM q 
is absolutely free of non-perturbative uncertainties. This is not the case in the old basis, 
which is therefore disfavoured. 

Thus we will use the new basis until very precise lattice values for all arising operators and the 
a s /mf,-corrections will be available. In 2007 also the Wilson coefficients of the corrections, 

i.e. have been determined^! Their size seems to be reasonably small. Since for some of 
the arising dimension eight operators even vacuum insertion approximation is not applicable, 
we do not include these terms in our numerics. 

3 HQE under attack 

Performing the HQE implicitly assumes that the sum over all possible exclusive final states, 
which is necessary to determine the total lifetime, is equal to the sum over all possible final state 
quarks. This assumption is called quark hadron duality, see e.gP^ or^Sl f or the application to 
decays of heavy hadrons. For total decay rates one obtains in the framework of the HQE 

r = r 0+ (A) 2 r 2 + (M 3 r 3+ (M 4 r 4 + .... (20) 

\m b J \m b J \m h J 

The leading term Tq describes the decay of the free b-quark and is therefore the same for all 
6-hadrons. It depends on the fifth power of the b-quark mass, so depending on what quark mass 
scheme one is using, one gets very different results. Therefore it is advantageous to look at ratios 
of lifetimes of different 6-hadrons. For different mesons this ratio starts with contributions of 
the order A 3 /mjj (for baryon vs. meson lifetime there are already corrections of the order of 
A 2 /m^), which look very similar to the contributions of T\ 2 0- Comparing lifetime predictions 
with measurements provides a test of the assumption of quark hadron duality. Although one 
should keep in mind that for T\ 2 one needs a stronger assumption, because less intermediate 
states (only the ones that are common to B q and B q ) contribute to the sum. To draw some 
definite conclusion on T\o also a direct determination of e.g. AT q is mandatory, which is now 
available for B s mesons^* 

In the last 20 years regularly some discrepancies between experiment and theoretical predictions 
assuming quark hadron duality have been arising, which always triggered a lot of theoretical 
interest. Below I will discuss the following topics: 

• In the mid 90ies the missing charm puzzle received a lot of attention, see e.g^ f OT a 
mini-review. The measured average number of charm-quarks per 6-decay turned out to be 
lower than theoretically predicted: n)? xp - < nf M . This discrepancy was also related to a 
mismatch of the experimental and theoretical values for the semi lepto nic br anching ratio. 
A possible solution was the violation of quark hadron duality, see e.g. I 65 | 66 | 



' A more detailed com parison of the theoretical determination of lifetime differences of different b-hadrons and 
of Y\ 2 is given e.g. in^3. 



• From the beginning of the 90ies the values for the lifetime of the A& baryon was measured 
to be much shorter, than the lifetime of the Bd mesons, while theory predicted them to 
be of similar size. This was also attributed to a failure of local duality in inclusive non 
leptonic heavy flavour decays, see e.g. 

• Before 2003 the values of T^ s /r^ d were smaller than 0.95. Here theory predicts almost 
exactly one. This would also hint to a violation of quark hadron duality or to new physics 
effects in the lifetimes. 

• In 2010^-^ and 2011^ the DO collaboration announced an unexpected large value of the 
di-muon asymmetry. It turned out that assuming new physics in Mf 2 is not sufficient to 
explain the large central value, one also needs large deviations of T\ 2 from its standard 
model value, either via a violation of quark hadron duality or via new physics effects, see 
e.g. EH. 

In the following we will show that the above puzzles disappeared or that the solution via a 
violation of quark hadron duality is not viable anymore. 



3. 1 Theory arguments in favour of the HQE 

In this section I will give some theoretical arguments in favour of the validity of the HQE. These 
arguments clearly represent no proof - for that an exact solution of QCD would have to be 
compared with the HQE, which is beyond the scope of the current paper. 

1. A pragmatic starting point is to calculate corrections within the HQE in all possible "di- 
rections" , in order to test the convergence of the expansion. Splitting up the contributions 
to Ar s in a leading term (no perturbative QCD corrections, vacuum insertion approxima- 
tion for the bag parameters and leading term in the HQE) and subleading corrections due 
to the use of lattice values for the bag parameters, due to perturbative QCD corrections 
and due to subleading HQE corrections one gets^J 

Ar s = Ar° (i + 5 Lattice + 5 QCD + <5 HQE ) 

= 0.142 ps" 1 (1-0.14 -0.06- 0.19) . (21) 

Explicit calculation shows that the corrections are at most 19%, which seems to be a 
reasonable value and indicates no breakdown of the theory. In particular the expansion 
parameter of the HQE is now determined to be 0.19 ~ 1/5. So no signal for a breakdown 
of the HQE is found by an explicit calculatior§. 

2. As mentioned above (see^lfor more details) the theoretical expressions for lifetime ratios 
of mesons resemble the expression for r \ 2 . 



(22) 



Moreover lifetimes are expected to be insensitive to new physics effects 0. Comparing 
theory and experiment for t(B + ) /r(-B^) one gets a very nice agreement^! Again, no sign 



9 This statement is of course much more solid than statements based on phase space arguments and simple 
power counting, see e.g. the discussion i m. 

h Sizable new physics effect in lifetimes correspond to sizable new physics effects in the dominant decay modes, 
which have not been seen experimentally. 



for a breakdown of the OPE. 

The precision of this comparison is however strongly limited by the corresponding hadronic 
parameters. The latest available lattice evaluation of the arising matrix elements stems 
from 2001^1 Here clearly an update would be very helpful. 

3. The inclusive approach to calculate Ar s described above, can also be compared to the 
exclusive approach, where one tries to estimate all exclusive decay channels that contribute 
to rf 2 . The seminal paper of Aleksan et. al from 1993^^ was recently update 

d E2] 

and the 

new results agree well with the predictions using the HQE. Again, an argument in favour 
of the validity of the HQE. 

We have summarized here several theoretical hints for the applicability of the HQE. These hints 
represent of course not a proof. In the next section we show how the above listed problems of 
the HQE predictions have been resolved experimentally. 



3.2 The missing charm puzzle 

The theoretical numbers for n c are more than ten years old^-^J so a reanalysis with updated 
input parameters seems to be desirable. Using a preliminary theoretical result with input 
parameters from 2012 and comparing them with the latest value from BaBar in 2006^^ one 
finds perfect agreement. 

2006BaBar _ -, on o+0.058 / O o\ 

n c — i.zuo_ 0055 . 

n s c M = 1.20 ±0.04. (24) 

So there is no hint for a violation of quark hadron duality from the charm counting anymore. 
Nevertheless it would be very desirable to have updated experimental numbers for n c from e.g. 
the full data sets of the B-factories, because this could shrink the allowed region for hypothetical 
new physics effects contributing to T\ 2 , see below. 

3.3 The A b -lifetime 

The experimental numbeP^ for the Ab-lifetime changed quite dramatically in the last ten years 

HFAG '03 r Ab = 1.229 ± 0.080 ps" 1 . (25) 
HFAG '12 r Ab = 1.425 ± 0.032 ps" 1 . (26) 

This is a shift of the central value of about 2.5o\ Part of this change is due to the fact that old 
measurements were only done with semi leptonic A& decays, while for newer measurements also 
non leptonic decays became available. The newer values can be easily explained by the HQE. 
e.g. while the old low values were very problematic. So again, the indications for a failure of 
the HQE have vanished. 

Looking a little closer, one sees however that there are still some issues that have to be settled. 
Experimentally there seems to be a discrepancy of a about 2 a between lifetime measurements 
of DO^and CDF^S]- both measure the non leptonic decay channel. To settle this issue we are 
eagerly waiting for new LHC results. 

But there are also some open theoretical issues, in particular the correct value of the non- 
perturbative matrix elements of the four-quark operators of the A^-baryon. Due to HQET only 
two different matrix elements (instead of four) arise. A widely used parameterization is 

■(A b \b Ll ,q L -q L rb L \A b )=:-^^r. (27) 



2m A * 48 



The second matrix element can also be expressed in terms of the parameter r, so in the end 
T A b / T B d — 1 is proportional to r. Unfortunately we do not have a lattice evaluation of r, only an 
exploratory study from 1999^^, whose numerical result should be taken with a lot of caution. 
In the literature values of r range from 0.2 (QCD sum rules^^) over 1.2 (the exploratory lattice 
study^^) up to 6.2 (QCD sum rules'^ - these results have met however a lot of scepticism in the 
sum rule community). Rosner related in 1996-^ the value of r to mass differences of 6-hadrons 
and obtained values between 0.9 and 1.8. Using his results and the recent PDG values^! of the 
6-hadron masses we get 



4 m|. - m| 6 
3 m B * — m B 



0.68 ±0.10, (28) 



which points towards smaller values of r and consequently to larger values of the A& lifetime. 
Here clearly new lattice determinations are mandatory to improve the theory prediction. 



3.4 The B s -lifetime 



For quite some time the B s lifetime was also considerably below the Bj lifetime. Just recently 
several new results for T s = 1/tb„ were obtained by performing the angular analysis in the 
decay B s — > ip4>. LHCbQSI CDfES and DO^ obtained from their fits: 



LHCb 
CDF 
DO 



t Bs = 1.520 ± 0.020 ps 

t Bs = 1.528 ± 0.021 ps 

1 A A q+0.038 ^ 
TB S = L44 3-0.035 P s • 



(29) 
(30) 

(31) 



Taking the most precise 

valuers] 

and the most recent standard model determination^^ we find 

an impressive agreement 



1.001 ±0.014 , 



T^SM 



0.996.. .1.000 . 



(32) 



This again confirms the framework of the HQE. For a more sophisticated comparison of course 
all lifetime measurements have to be combined properly - this was done also recently by HFAG 
(for this average the most recent LHCb numbeff^ is not yet included, but the preceding 
number fromLLZ') : 

— = 0.984 ±0.011 . (33) 



This new average agrees within about 1<t with the very precise theoretical prediction and the 
HQE seems to work with a high accuracy. 

The new lifetime values for the B s meson can also be used to update the theory predictions 



for the effective lifetimes. Previously 1 ^ t Bs = L477Io;o22 

ps was used. Replacing this number 
by the most recent LHCb measurement^ one finds the new SM model predictions for the 
effective lifetimes and the flavour-specific lifetime listed below. This can be compared with 
the experimental numbers for r EfI (K + K~) from LHCb^% for T Efr (V>/ ) from CDF^H and the 
average of the flavour specific lifetime from hfagISI 





Exp. 


SM-old 


SM-new 


T* n (K + K-) in ps 


1.468 ±0.046 


1.390 ±0.032 


1.43 ±0.03 


T m (ipfo) in ps 


1.70 ±0.12 


1.582 ±0.036 


1.63 ±0.03 


T bti in ps 


1.463 ±0.032 




1.54 ±0.03 



For the SM predictions of the effective lifetimes the dominant sources of the error are the theory 
value of Ar s (about ±0.02) and the experimental values for tb s - The effective lifetimes agree 
within la, while the experimental value for the flavour-specific lifetime is about 2.4<r below the 
theoretical central value. Here also new, more precise data will be very desirable. 

3. 5 First measurement of Ar s 

In Moriond LHCb presented the first measurement (> 5<r) of a non-zero value of the decay rate 
difference in the neutral B s system^S! 

Ar s = (0.116 ± 0.019) ps" 1 . (34) 

Naively this corresponds to a statistical significance of 6.1a. AT S was obtained from the angular 
analysiP^l in the decay B s — > J/%[j<fi. CDP^and also performed similar analyses to get: 

DO 81b" 1 Ar s = (0.163 ± 0.065) ps" 1 , , . 

CDF 9.61b" 1 AT S = (0.068 ± 0.026 ± 0.007) ps" 1 . ^ ' 



The new HFAG average^^ reads (for this average the most recent LHCb number^ is not yet 
included, but the preceding number fromLLTJ) : 

Ar s = (0.100 ± 0.013) ps" 1 . (36) 

There is also a long history of theoretical predictions for Ar s . The first estimates stem from 
1977EB where AT d /T d = 1/6 was found (for B d mesons!), which is very close to the current value 
for Ar s /r s . As can be seen from Eq . (|21j) NLO-QCD corrections, subleading HQE corrections 
and the impact of lattice parameters for Ar s are ample. Therefore the numerical values of old 
calculations have to be taken with a pinch of salt, even if some of them are by accident close to 
the current experimental number. 

The subleading l/m^-corrections were determined in 1996"" j The dominant NLO-QCD correc- 
tions were calculated in 1998^, subdominant ones in 

2003 [13114J In 2006^21 

more theoretical 

improvements were included - now also the relatively well-known bag parameter Bb s gives the 
dominant contribution. Before that the dominant contribution came from the bag parameter 
Bs, which describes the matri x element s of the operator Qg, see. Eq, (|16p . First lattice estimates 
of B$ appeared from 2000 on l 88 | 89 [90]g o ^ before 2000 in principle no reliable estimates of Ar s 
could be given at all, because of missing inf ormati on and before 2006 estimates were affected by 
large uncertainties. More recent evaluations ^^ ^ read § 

2006 : Ar s = (0.096 ± 0.036) ps" 1 , (37) 
2011 : Ar s = (0.087 ± 0.021) ps -1 . (38) 

We will stick to Eq. (]38p as our theory prediction for AT S . Comparing experiment with theory 
we again find a perfect agreement: 

Ar? xp - 0.100 ±0.013 



Ar? M 0.087 ±0.021 



1.15 ±0.32. (39) 



This result proofs experimentally that the HQE works also for the calculation of T? 2 • The issue 



12- 



of the proper calculation of T\ 2 gave raise to numerous discussions in the literature, see e.g. 
and references and citations of this paper. This issue is now settled. 
Now several comments to Eq. (|39|) are appropriate: 



"The Rome Group will soon present^^ a numerical update, which yields Ar s /r s = 0.149 ± 0.015. This 
corresponds to AT S = 0.098 ± 0.010 ps -1 and is in perfect agreement with the measurement. It also agrees with 
the result from^4 



Eq. (|39p shows that the HQE works for rf 2 with an accuracy of 15% ± 32%. Therefore the 
question is not anymore, whether the HQE is appropriate for rf 2 , but how precise is the 
HQE. Does it work to an accuracy of 30%? Or maybe even to an accuracy of 10%? This 
can be answered by further reducing the theoretical and the experimental error. 

The uncertainty in Eq. (|39p is dominated by the theoretical error, but an improvement of 
the experimental error will also be helpful. Here new data from LHCb and the two planned 
Super-B-factories ^^ l ^ will be very important. 

The theory prediction for T\ 2 has again a crucial dependence on non-perturbative pa- 
rameters. In particular it is also proportional to f B . For the theory prediction we have 
used again fg g = 231 ± 15 MeV. As in t he case of Mf 2 we also show the results, if new 
evaluations of the decay constant are used^ES 



Reference 


fs s 


Ar s 


HPQCD (1110.4510) 
Fermilab/MILC (1112.3051) 


225 ± 4 MeV 
242 ± 10 MeV 


0.083 ± 0.017 ps" 1 
0.096 ±0.021 ps" 1 



(40) 



In contrast to M\ 2 now the error is not reduced so strongly, if one uses smaller errors for 
the decays constant. For the same reasons as discussed in relation to AM S we will use the 
prediction from Eq. ([38|) . 



To see, what is necessary to increase the theoretical accuracy, we compare the error budget 
from the theory prediction of AT S in 2006^21 with the one from 201 ffH 



A pbM 


2011 


2006 


Central Value 


0.087 ps" 1 


0.096 ps" 1 




17.2% 


15.7% 


S(f Bs ) 


13.2% 


33.4% 


SQj) 


7.8% 


13.7% 


S(Bs,b s ) 


4.8% 


3.1% 




3.4% 


3.0% 


S(V cb ) 


3.4% 


4.9% 


S(B Bs ) 


2.7% 


6.6% 










24.5% 


40.5% 



Currently the dominant theoretical uncertainty stems from the power suppressed operator 
i?2- These matrix elements have currently only been estimated with vacuum insertion 
approximation. Nu meric al subleading effects within the QCD sum rule approach have 
been determined inEHES_ Here a more precise determination of Br 2 is mandatory in order 
to improve the theoretical accuracy of T^ 2 - 

For a long time the dominant error source of the theory prediction of Ar s was Jb s , currently 
the dependence on f B is the second largest error. Thus one might want to get rid off the 
dependence on fs a by considering the ratio Ar s /AM S - here one implicitly also assumes 
that no new physics acts in AM S . In this ratio the overall factor f g cancels, as well as 
the bag parameter Bb„ of the operator Q, given in Eq.flSJ): 



Ar s 

AM, 



B' 

46.2 + 10.6-f 
B 



10" 

0.0050 ±0.0010 . 



13.2- 



B 



R2 



B 



2.5 



B Ro 
B 



1.2 



Br 
B 



(41) 



The numerical dominant term is now completely free of hadronic uncertainties, the remain- 
ing dependence on hadronic parameters is always via ratios of bag parameters. Comparing 
this theory predictionHD with the experimental number presented by HFAG^31 one g e t s 

/ Ar s \ Exp I ( Ar, ^ SM 



{aWJ = i- 1 ^ - 27 - ( 42 ) 



Once more an impressive agreement with the HQE prediction and the assumption that no 
new physics acts in AM S or AT S . 

3.6 Ar s from Br(B s -> D { * )+ + D^~)9 

In the literature one finds regularly the following relation 

Ar?p 



r 



s 



2Br (B s -> D^ )+ + £>M~) (43) 



and one might of course be tempted to use Eq. (|43p . to determine ArJ? p = |rfg| directly from 
"simple" branching ratio measurements. This equation was derived in 1993 to show the 
equivalence of the exclusive and inclusive approach for calculating Ar s . It holds in the limit 
m c — > oo; rrif) — 2m c —> 0;N C — > oo. This limit corresponds to neglecting the 3-body final state 
contributions to Tf 2 - In 1993 Aleksan et al. estimated^ 

Ar 

oc 0(0.15) , (44) 

1 s 

which is amazingly close to the current experimental value ^31. This analysis was redone 
recently^. With up to date values of the input parameters one gets: the 2-body final states 
contribute 0.100 ± 0.030 to Ar s /T s . This number is the updated value of Eq.(|44|). It is slightly 
below the HQE prediction in Eq. (|38p . Moreover it turned out that the 3-body final state con- 
tributions - which are neglected in the derivation of Eq. (l43l) - are quite sizable, they contribute 
about 0.06... 0.08 to Ar s /r s . Since this is comparable to the 2-body final states, the above 
approximation turns out to be a bad one. We repeat here the suggestion fronts 
We strongly discourage from the inclusion of Br(Bs — > Z)M+ + Z)M~) in averages with AT S 
determined from clean methods. 



3. 7 HQE at work: semi leptonic CP asymmetries and mixing phases 

In the previous sections we have shown that the HQE describes inclusive decay of b-hadrons very 
accurate. Even the prediction of rf 2 agrees with an accuracy of about 30% with the experiment. 
To make more precise statements, some non-perturbative progress is mandatory. 
Because of this success we use now the HQE to predict^ also quantities that might be very 
sensitive to new physics effects like semi leptonic CP asymmetries and mixing phases. 

a s fs = (1.9 ± 0.3) • 10" 5 , <p s = 0.22° ± 0.06° , (45) 
a^ = -(4.1±0.6)-10~ 4 , 4> d = -4.3° ± 1.4° , (46) 
A b sl = 0.406afj + 0.594af; = (-2.3 ± 0.4) • 10 -4 . (47) 

The semi leptonic CP asymmetries turn out to be extremely small, as well as the mixing phase 
in the 5 s -system. Therefore any measurement of a sizable value would be a clear indication of 
new physics. The so-called di-muon asymmetry A b sl is a linear combination of the semi- leptonic 
asymmetries in the Bj and B s systems. As mentioned in the the introduction there was a 



considerable excitement in the past, because we had several hints for huge new physics effects 
in 5-mixing. Before the Lepton-Photon conference 2011 a frp^of flavour observables resulted in 
a huge deviation of the B s mixing phase from the tiny SM prediction: 



-SM 



s 



+ <t>f , (48) 
cf)f = -51.6° ± 12° , (49) 

where (f)f~ denotes the de viatio n from the SM value of the l^-mixing phase. Moreover the 
DO collaboration measureJP^^ a large deviation of the di-muon asymmetry from the theory 
expectation: 

A b sl = -(7.87 ± 1.72 ±0.93) • 10" 3 . (50) 

This number is by a factor of 34 larger than the SM expectation, the statistical significance of 
the deviation is 3.9a. The di-muon asymmetry currently still has its large value, while the hints 
for large deviations of <f> s from zero have disappeared. A fit of the most recent data, e.gff^ from 
B s — > J/ip<j) gives 

<fi = -0.1° ±6.1°. (51) 

This number is perfectly consistent with the SM, while still some ample deviations are allowed. 
Current experimental bounds on the individual semi leptonic CP asymmetries are given by 
HFA(P] 



a s fs = (-1050 ± 640) • 10 -5 , (52) 
a d fs = - (33 ± 33) • 10" 4 . (53) 

These remaining deviations will be investigated systematically below. 

3.8 The decay rate difference in the B^-system 

The SM predicts^^ also a tiny width difference in the B^ system. 



(4.2 ± 0.8) • It)- 6 . (54) 



Fa 

This can be compared to the current most precise bounds from Belle. 



(-17 ±21) • 10~ 3 . (55) 



Ar s is governed by the CKM- leading tree-level decay b — > ccs. So any sizable new physics effect 
to this decays would also affect strongly many other flavour observables. Ar^ is governed by 
Cabibbo-suppressed decays, hence the bounds on new physics effects in Ar^ are weaker and the 
real value of A^ might deviate from its almost vanishing SM value. This interesting null-test of 
the SM was discussed recently in^l 

4 New Physics in S-mixing 

New physics in i?-mixing can be paramterized model independently in the following way 

T\ 2 = T^ SM , M? 2 = M?£ M -A q , A q = | A, | e< . (56) 

Here we assume that new physics acts only in Mf 2 , while T\ 2 is given by the SM prediction. 
Although this is expected to hold within the hadronic uncertainties (up to 25%), we will also 
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Figure 2: Bounds on the complex Aq-plane from AM, (red), AT q / AM q (yellow), a q sl (turquoise) and direct 

determinations of <j> q (black rays). 



investigate new physics effects in T\ 2 in the next section. 

Using the parameterization of Eq. (|56p one can express the observables in the mixing system in 
terms of the SM predictions for M\ 2 and T\ 2 and in terms of the complex parameter A q : 

|A,| , (57) 
>f + <) , (58) 



AM q 


= 2|M x f M |- 


AT q 


= 2|r| 2 | • cos 


AT q 


|r? 2 | « 


AMg 


\M?f M \ 




|r? 2 | s 

\Mg M \ 



(< 



(OS + 



IM 



sin (</> 9 SM + 



|A 9 | 



(59) 
(60) 



Noticing sin(^g M ) ~ 1/240, one sees that new physics easily could enhance the semi leptonic 
asymmetries in the -B s -system by a factor of almost 250, if a sizable phase is present in Mf 2 . 
The relations of Eq. ([57]) - Eq. ([60j) can be used to obtain bounds on the complex Ag-plane. As 
inputs one needs the SM predictions for M^ 2 and r^ 2 , as well as the experimental numbers for 
AMg, Ar 9 , a q sl and (j> q = <^ M + </> A . For the hypothetical case |A g | = 0.9 and 4>f = -7r/4 one 
gets the bounds of Fig. ([2]) in the complex A g -plane. Eq. (i57"|) gives a bound on the modulus 
of A q , this corresponds to a circle in the complex A q -plane (red region in Fig. ([2])). Eq. (|58p 
gives a direct bound on the phase, this corresponds to the black rays in Fig. ([2]). Because of 
the relatively large theory uncertainties in r| 2 , it is however much more advantageous to use a 
direct determination of (p q , e.g. for <f> s from the angular analysis of decays like B s — > J /ipcj) or 
B s — > J/i/jflf^. Eq. (j59j) and Eq. (|60l) constrain simultanously the modulus of A q and the phase 
A . AT q / AM q gives the yellow region in Fig. ([2]), while a q sl gives the turquoise one. 
Combining all data till march 2012 and allowing only for new physics in Bj- and i? s -mixing 
(i.e. neglecting e.g. new physics penguin contributions), one gets the following bounds on the 
complex A^- and A s -planeP^ (update o ftH). 



excluded area has CL > 0.6B 
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New Physics in - B d mixing 



excluded area has CL > 0.6B 
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In 20l|Q new physics in l?-mixing could very well accommodate the different deviations from 

the SM expectations, seen at that time. This is not the case anymore in 201^^. There is now 

a tension between the direct determination of <p s and the di-muon asymmetry. 

In the .B^-system, new physics in M\2,d can resolve the discrepancy between B —> tv and direct 

determinations of sin 2(5. In the -B s -system everything looks SM-like although still sizable values 

for (j)^ are possible. Just recently a second (symmetric) solution in the complex A s -plane was 

excludedP^. We also would like to note that in^Sl n0 tension is found for ex ■ 

To improve further the bounds on the complex A„-planes, more precise data are necessary. 



5 New Physics in Ti2 

The theory expression for the di-muon asymmetry can be written in the following way 
A sl = (0.594 ± 0.022) (5.4 ± 1.0) • lO' 3 ^™ + ^) 



■ sin 



|Ad| 

hSM 



+(0.406 ± 0.022)(5.0 ± 1.1) • 10~ 3 ^^- ' ^ s ' . (61) 

\^s\ 

Since A s and are bounded from measurements of the mass differences to be close to one and 
the sine can be at most one, there exists a theoretical upper limit for the di-muon asymmetry. 
We use here the fit values of A q from^Slto obtain the following upper bounds: 



A s i < 



-1.7 


lO" 3 : 


la for 


|A q 


-2.8 


10^3 : 


3a for 


|A 9 


-7.5 


10-3 : 


3a for 


|A 9 



la for <fiq, 

3a for <j)f, (62) 
set sine to 1. 



For the first number the four parameters of A„ (q=s,d) have been chosen to take the value, 
which gives the largest di-muon asymmetry, within the allowed 1 a range of the fit in^Sl ) f or 
the second number, the 3 a range has been chosen, while for the third number the sine has been 
set to one by hand. The last number is purely hypothetical, because such a large value of the 
mixing phase is in contrast to experimental investigations of e.g. B s — > J/ip(j)%. The above 



J This also holds, if one takes into account large new physics penguin contributions to the decay b —¥ ccs, which 
could lead to a certain extent to a cancellation between the penguin phase and 4>q ' ■ See the discussion in the next 
section. 



bounds have to be compared with the experimental measurement 



A%° = (-7.8 ±2.0) • 1(T 3 . (63) 

One fin ds th at the central value of the di-muon asymmetry is larger than theoretically possible, 
see also^ ^. This discrepancy triggered a lot of interest in the literature (see the list of citations 
f or | 6|7|8 |j rp^g theoretical upper limit stems from the fact that we did not allow ample (i.e. 
larger than the hadronic uncertainties, which are of the order of 25 %) deviations of T\ 2 from its 
SM value. One might of course be tempted to give up this assumption. If there is no new physics 
acting in M\ 2 one would need an enhancement factor for Y\ 2 of 34 to describe the central value 
of the di-muon asymmetry, if new physics acts in Mf 2 one still needs an enhancement factor of 
about 4.6 (1 orange of Eq.([62D) or 2.8 (3 a-range of Eq.flM}). 

As discussed above violations of quark hadron duality in the range of 280% — 3400% are clearly 
ruled out now, as experiment and theory agree for Ar s to an accuracy of about 30%. 
Next one might think of new physics effects in T\ 2 . One can modify Eq. (j56p to allow for new 
physics in T\ 2 : 

T\ 2 = Tli M ■ A q , M\ 2 = Mti M ■ A q , 
A q = \A q \e^ q , A q = \A q \e<. (64) 

Just looking at AM q , AT q and a , one can of course fit any numerical value of these three 
observables with the three parameters | A„ | , \5„\ and A ^ — ,,s 
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|A,| 



(65) 



,SM + A _ A g (66) 





(Vf + < - ^) . (67) 

So the situation AM q = AM™, Ar g = AT s q M and a% = 34a^ SM can be described by |A 9 | = 1, 

| A, | • cos (d% M + (j) q - (f>~ Aq ^ = 1 and |A 9 |/|A 9 | • sin [cjyf 1 + </)f - (jfiA =1, which is equivalent 

to |A| = 1, |A| = 34.0147 and cos (<^ M + (j) q - </A) = 0.0293991. In this spirit bounds on new 

physics contributions to Tf 2 and Tf 2 were obtained in^l, see Fig.©. It turned out to be useful 
to introduce a new parameter 5 q : 



1 12 



*. = / '.'sm x ■ (68) 

Ar / / Ar SM \ 
Re < 5 «> = aW,/{a^) • < 69 > 

toft) = -<,/{^m) ■ (TO 

The experimental constraint on Ar s can be fulfilled by Re j s ~ 1 and the semi leptonic asym- 
metries only affect Im 5 S . The real part of 5^ is almost unconstrained: Re 5^ ~ —4.3 ± 5.4 (see 
Section T3.8p . while the the semi leptonic asymmetries again affect only the imaginary part of 84. 
Y\ 2 is dominated by the Cabibbo- favoured decay b — > ccs. Thus any sizable contribution to 
Tf 2 will affect many other observables dramatically. As explained in^^l these effects on other 
observables have to be taken into account, when using the results of FigQ! Any new physics 




Figure 3: Bounds on the parameter lm(S q ), if only AM,, Ar, and a q al are taken into account. Large values of 
Im(<5 s ) violate however the perfect agreement of experiment and theory for quantities like rs s /r_B d ,... (see text). 



contribution to Tf 2 corresponds to a transition bs —> X, where X can be one or more (SM or new 
physics) particles and Mx < Mb s - The new operator (bsX) gives now rise to many different 
new contributions. Via diagrams like (if X corresponds to two particles) 




the (6sX)-operator contributes to the leading term T in the HQE, see. Eq, (f20]) . This affects 
the total decay rate T s and therefore quantities like the semi leptonic branching ratio, or the 
branching ratio of a 6-quark into two, one or zero charm quarks (see the discussion of the missing 
charm puzzle in Section( j3.2p ). Via diagrams like 



b b 




the (6sA)-operator contributes to the third term T3 in the HQE, see. Eq, (|20p . This affects 
the lifetime ratio t(B s )/t(B c [), but also the decay rate difference AT S and the semi leptonic 
CP-asymmetrief|. It will also give rise to new contributions to Mf 2 . Via mixing such a new 
operator can also modify well constrained quantities like the rare decay b — > s"f. 
All in all, any new physics contribution to Tf 2 will face severe constraints from many well 
measured quantities. Hence it seems impossible that ample (i.e. larger than the hadronic 
uncertainties, which are of the order of 25%) new physics effects are present in rf 2 . Unfortunately 
this statement cannot be proven model independently. At least the operator structure of the 
new contribution has to be specified, so that one can also calculate the new contribution to e.g. 
To and l^V Such an analysis'^ has been done for the most promising (because least constrained) 
candidate, X = t + + r~ and it turned out that new physics contributions to rf 2 have to be 
smaller than 30% or 40%, depending on the Dirac-structure of the 6srr-operator. This clearly 
excludes the possibility of having new physics contributions to Tf 2 in the range of 280% — 3400%. 
More precise data on rs a /rB d , B s i, n c , Br(B s — > t + t~),... will further shrink the bound on new 
physics effects in Tf 2 (or give hints for small new physics effects). 

Finally there is also the possibility that DO simply saw a 2.5 a upward fluctuation of the result 
for the di-muon asymmetry and that the actual value is below -2.8 per mille, which still would 
be a large deviation from the SM. To get a final clue about that, independent measurements of 
the semi leptonic asymmetries are urgently needed. 



6 Penguin pollution 

In the last section I will discuss the question to what extent does the value of the di-muon 
asymmetry contradict the result of analyses of B s — > J/ip(ft from LHCb and Tevatron. The 
angular analysis of the decay B s —> J/ip(j) at CDF, DO and LHCb gives among others the 
quantity S^. If one neglects certain penguin contributions, then is given in the SM simply 

fc Strictly speaking the contributions to and Fs are of course not equal, see e.g. the discussion in^l 



by CKM-elements. Its SM expectation readiSP^ 



c-SM 



0.036 ± 0.002 . 



(71) 



In the presence of new physics will be modified by the new mixing phase 4>f~, which also 
changes the values of the semi leptonic CP-asymmetries 



Q SM 

D 1p(f> 



sin 



(2/9.-0 



0.00 ±0.11 



(72) 



As already mentioned above, in deriving e.g. Eq. (l72p certain pen guin contributions have been 
neglected. Since the current experimental bound on (taken frorrA^j) has already a small error, 
this approximation might not be justified any more. It was already pointed out in^^H that a 
measurement of a relatively small value of the mixing phase requires a much more stringent 
control of the penguin contributions. Including penguins one gets the general relation 



cSM 

D lpcj> 



sin 



(2&-<A 



A _ ^Peng,SM _ ^Peng,NP 



(73) 



Looking at Eq. (|73p one sees that there is the possibility that a non- vanishing value of the new 
mixing phase <f>^ is compensated by SM and new physics penguins. Now it depends on the 
possible size of the penguins, to what extent such a cancellation can occur. 
To exactly determine the size of the penguin contributions one would have to solve all the non- 
perturbative dynamics related to this decay, which is clearly not possible. For a rough numerical 
estimate we start with the amplitude 



A(B S ^ f) = (B s \n ef f\f) 
with the effective SM hamiltonian for b — > ccs transitions 



H 



eff- 



Gf 
V2 



Xu {CiQl + C 2 Q U 2 ) + A c {C x Ql + C 2 Q C 2 ) + X t £ dQi 

i=3 



(74) 



(75) 



The CKM structure is given by \ x := V* b V xs ; the decay b — > ccs proceeds via the current-current 
operators Qi,Q 2 and the QCD penguin operators Q^,...,Qq. C\,...,Cq are the corresponding 
Wilson coefficients. When the current-current operators Qi,Q 2 are inserted in a penguin di- 
agram in the effective theory, they also contribute to b — >■ ccs. Electro- weak penguins are 
neglected. 

Therefore we have the following structure of the amplitude A (B s — > f) 



A 



Gf 
V2 



u\T 
i I 



i=l,2 



i=l,2 



i=3 



(76) 



(Q} T denotes the tree-level insertion of the local operator Q, (Q) p denotes the insertion of the 
operator Q in a penguin diagram. Using unitarity At = — X u — X c we can rewrite the amplitude 
in the following form. 



A — —L\ 



E c i(Q 



c\T+P 



1=1,2 



i=3 



A 



u\T 

i , 



,1=1,2 



i=3 



.(77) 



The first term gives rise to f3 s in the CP asymmetry. The second term (proportional to X u ), 
corresponds to the SM penguin pollution. To get an idea about the possible size of the second 
term we notice 1 ^ that Re (A u /A c ) < 0.01 and Im (A u /A c ) < 0.02. Moreover the penguin 
Wilson coefficients IC3... q\ are typically smaller than 0.04, therefore one can neglect them in 



comparison to the Wilson coefficients C\,Gi ~ 1. For inclusive decays it turned out that 
(Q) p < 0.05(Q) T , which would result in a penguin pollution smaller than 0.0005 + O.OOlOi, i.e. 
smaller than one per mille compared to the contribution of the golden plate mode and thus 
completely negligible, even if one takes into account the current experimental precision. Of 
course this estimate is very naive and it might change considerably if the matrix elements of the 
exclusive final states are taken into account, but also a non-perturbative enhancement of this 
estimate by a factor of 10 would not change the conclusion. 

Here clearly more theoretical work has to be done to quantify the size of penguin contributions, 
Experimental strategies to dete rmine the penguin contributions to B s —> ipcf) and B s — > tpfo 
have been discussed in e.g. QUE anc j ^ ne second reference of^l An investigation of the decay 
B s — > ^Ke; at t he LHCb upgrade can also gain important insights on the size of penguin 
contributions ^21. Nevertheless it seems that the penguin pollution in the SM is tiny, except 
some unknown non-perturbative would arise. 

In principle new physics penguins might be much larger, but it is hard to imagine that a phase 
(f)f~ of order one can be compensated by penguins, without violation other constraints. Therefore 
the angular analysis of B s —> tpcp seems to be in a slight contraction to the central value of the 
di-muon asymmetry. 

Finally we will show that even small penguin contributions have an observable effect and that 
they also can be analysed quite model independent The mass difference, the decay rate 
difference, the semi leptonic CP asymmetries and S^<p are related via 
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All SM phases are small: ^ M = 0.22° ± 0.06° and -2f3 s = (2.1 ± 0.1)°. Nevertheless they are 
non-neg ligiblPl This can be seen by drawing 5 in dependence of 4>^ for different values of the 
penguin contributions, JP en g< SM + <5P en s> NP = 0°, 2°, 5°, 10°. 
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The curve closest to S = 1 corresponds to £P en g> SM + <5P en s> NP = 0°, the next one to 2° and so 
on. Because of this huge sensitivity to penguin effects Eq, (|78p can be used to gain some insight 
in the size of penguin contributions. 



7 Conclusion 

In Moriond 2012 the first measurement of AT S was presented from LHCb. Combining the LHCb 
data with the corresponding ones from CDF and DO and comparing their average with theory 
predictions one finds an impressive agreement: 

ATg 1 = 0.100 ± 0.013 = 

Arf M 0.087 ±0.021 v 1 

This represents the first experimental proof that the HQE can also be applied to the calculation 
of Tf 2 , with an uncertainty of about 30%, which is an important theoretical insight, because T\ 2 
is expected to be most sensitive to violations of quark hadron duality. To some extent this is 
quite amazing, because the energy release in the dominant decays, contributing to Ar s is not so 
large: mB s — 2m£> s ~ 1.428 GeV. Currently the experimental and the theoretical errors for Ar s 
are of similar size. To improve the theoretical accuracy, matrix elements of higher dimensional 
operators have to be determined non-perturbatively. 

Total decay rates, i.e. lifetimes, are expected to be less sensitive to violations of quark-hadron 
duality. There are now several precise numbers for t(B s ) available from CDF and LHCb, which 
also agree perfectly with theory predictions. 

Tr Ex P Tr SM 

= 1.001 ±0.014 , = 0.996.. .1.000 . (81) 

T B d TB d 

This is again a perfect confirmation of the applicability of HQE. In the case of the lifetime 
ratios tr+ l T B d and T\ b /rB d no deviation is currently existing, the precision of any comparison is 
however strongly limited by our poor knowledge of the corresponding hadronic parameters and 
also by the uncertain experimental value of T\ b . 

In accordance with this perfect agreement between experiment and theory, the new LHC and 
Tevatron data show no further hints for new physics in B s -mixmg. A model independent fit 
of all flavour data is consistent with no new physics in B s -mixing (although there is still some 
room for sizable deviations from the SM expectations) and some small deviations in .B^-mixing. 
To investigate these issues further a better control over penguin contributions is mandatory. 
The large central value of the di-muon asymmetry is still an unsolved problem, because it can 
not be explained by new physics contributions to Mf 2 alone, instead one needs in addition an 
enhancement of Tf 2 compared to its SM value by a factor of 2.8 up to 34. In the text is was 
discussed, however that rf 2 can not deviate more than about 30% from its SM value (neither 
due to new physics nor due to violations of quark hadron duality). On the other hand one 
should keep in mind that the central value measured by DO is only about 2.5 a above the 
theoretical limit. To settle this problem clearly an independent measurement of semi leptonic 
CP asymmetries might shed light into the dark. 

The success of the HQE for predicting inclusive processes in the 6-system, raises of course again 
the question of the applicability of the HQE to the charm system. Which seems now at first 
sight to be more promising, because for D-mixing the energy release is not so different from 
Ar s : mo — 2m,K ~ 0.9 GeV and mo — Im^ ~ 1.6 GeV. Of cours e the situation is much more 
subtle because of the huge GIM cancellations in L>-mixing, see e.g. I 1Q*H 
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